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We propose a scheme for a quantum thermal machine made by atoms interacting with a single 
non-equilibrium electromagnetic field. The field is produced by a simple configuration of macroscopic 
objects held at thermal equilibrium at different temperatures. We show that these machines can 
deliver all thermodynamic tasks (cooling, heating and population inversion), and this by establishing 
quantum coherence with the body on which they act. Remarkably, this system allows to reach 
efficiencies at maximum power very close to the Carnot limit, much more than in existing models. 

Our findings offer a new paradigm for efficient quantum energy flux management, and can be relevant 
for both experimental and technological purposes. 


I. INTRODUCTION 


Recent years have seen an uprising interest in ther¬ 
modynamics at atomic scale m due to the latest- 
generation manipulation of few, if not single, elementary 
quantum systems [M]. In particular, out-of-equilibrium 
thermodynamics of quantum systems represents one of 
the most active research areas in the field |MQ]. In this 
context, triggered by vast technological outcomes mm, 
the concept of quantum absorption thermal machine [T5] 
has been reintroduced [TH - THH . These machines are par¬ 
ticularly convenient since they function without external 
work, extracting heat from thermal reservoirs through 
single atomic transitions to provide thermodynamic tasks 
(e.g., refrigeration). 

Nonetheless, fundamental issues remain unsolved. A 
first one is the connection a single atomic transition to 
given thermal reservoirs, posing serious obstacles to prac¬ 
tical realizations of such machines. A second, more theo¬ 
retical issue concerns the role of quantumness. Indeed, in 
typical models quantum features are not required mm, 
and only recently the advantages of quantum properties 
in thermal reservoirs have been pointed out |16j . The 
role of quantum features in the machines itself is debated 
mm, so that the advantages of quantum machines over 
standard ones remains partially unclear. 

In this paper we address both of these open prob¬ 
lems by introducing a new quantum thermal machine 
setting, based on an out-of-thermal-equilibrium (OTE) 
electromagnetic bath naturally (i) coupling to each single 
atomic transitions, and (ii) creating quantum features in 
the machine. The field is produced by macroscopic ob¬ 
jects, and acts on each atomic transition as a different 
thermal bath at an effective temperature, hence provid¬ 
ing all the elements needed for quantum absorption tasks. 

This paper is structured as follows: the physical system 
is introduced in Section [TTJ along with the master equa¬ 
tion governing its dynamics, while Section m is devoted 
to the introduction of thermodynamic quantities char¬ 
acterising the heat exchanges happening between atoms 
and field. In Section [TV] the first part of the results is 
given, concerning the action of the machines, the dif- 



(a) 



FIG. 1. (a) A slab of thickness 5 at temperature Ts is placed 
in the blackbody radiation of some walls at temperature TV- 
Two atoms are placed in the resulting OTE electromagnetic 
field, at a distance z from the slab and at a distance r from 
each other, (b) Stationary heat fluxes between the OTE en¬ 
vironment E and each atomic transition. The OTE held also 
mediates for an effective atomic interaction producing an en¬ 
ergy exchange Q r between resonant atomic transitions. Each 
flux contribution corresponds to a term in the atomic master 
equation |l|. 


ferent tasks it can produce and their intrinsic quantum 
origin. The second part of the results of this work about 
the machine efficiency and its Carnot limit is given in 
Section |V[ Finally, remarks and conclusions are drawn in 
Section iVTl 
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II. PHYSICAL SYSTEM 

The setup of this paper is schematically depicted in 
Fig. 0 where a slab of thickness 6 at temperature Tg 
is placed in the blackbody radiation emitted by some 
walls at temperature T\y 7 ^ Tg. The total electromag¬ 
netic field embedding the space between the slab and the 
walls is therefore given by the sum of four contribution: 
the direct blackbody radiation of the walls, the radiation 
emitted naturally by the slab and the walls’ radiation 
after being either reflected or transmitted by the slab. 
Such an OTE field has been studied in the context of 
Casimir-Lifshitz force and heat transfer [5DH2S1 1 where its 
properties have been characterised in terms of the field 
correlators through a scattering matrix approach. The 
slab and the walls, macroscopic objects, are here the only 
ones directly connected to thermal baths. In addition, a 
three-level atom M (machine) and a two-level atom B 
(target body) are placed at the same distance 0 from the 
surface of the slab and spatially separated by a distance 
r. The atomic open system involves then four transitions: 
the body transition labeled as B and the three machine 
transitions labeled as 1,2,3. Transition 1 connects the 
two l owest-lying energy eigenstates (red transition in Fig. 
|l(b)[ ) and transition 2 the two highest ones (green tran¬ 
sition in Fig. |l(b)[ ). The OTE field interacts with them 
through the Hamiltonian Hj = — JT d* -E(Rj), where d, 
is the dipole moment of the i-th transition of the atomic 
system and E(R,) is the electromagnetic field at its po¬ 
sition R,. The total Hamiltonian of the system is 

Htot = Hm + Hb + Hfield + H[. (1) 

where iFfieid is the Hamiltonian of the OTE field. In 
the following we will not need the explicit expression of 
iffieid since only the field correlations will enter the mas¬ 
ter equation describing the dynamics of the atoms. The 
free atomic Hamiltonians Hm and Hb have expressions 

= (w_B + AS , (w_B))u- t cr = WgerV, (2) 

H m = (wi + AS(ui) + S~ (w 2 ) - S _ (w 3 ))/4ki 

+ (w 3 + A S(u 3 ) + S + (u} 2 ) - S~(u i))k 3 k 3 

= wi«}ki + w 3 k^k 3 , (3) 

being a (cr) the lowering (raising) operator of the body 
B and K n . K' n (n = 1, 2,3) the lowering and raising oper¬ 
ators of M. S ± (u}i) here represents a shift of the energy 
of each level in the i— th transition due to the local inter¬ 
action with the field and A S(w) = S + (uj) — S~{w). In 
the third equality the renornralised transition frequen¬ 
cies u>i for M and B have been introduced to account 
for the effects of the shifts S ± (iOi). Throughout this 
work we will always assume the physical consequences of 
such frequency renormalisation to be negligible, such that 
uii = LOi Vi. This assumption has been fully confirmed by 
extended numerical simulations, having always detected 
the relative error introduced by neglecting these shifts to 
be less than 1%. 


It is worth stressing here that, differently from previous 
works on atomic-scale thermal machines [HnTB; ITS . T9] , 
each atomic transition interacts here with the same elec¬ 
tromagnetic field, which embeds all the space where the 
atoms are placed. As we will show in what follows, there 
is then no need to conceive different environments, each 
interacting with a single atomic transition: a single non¬ 
equilibrium electromagnetic field is here able to produce 
all the physics needed for quantum thermodynamic tasks. 


A. The master equation 

In [24] the master equation (ME) for two emitters in 
such a field has been derived under the Markovian limit 
as 

f P\ + D b (p) + D n{p) + D d (p), (4) 

n =1 

where H T = H M + H B + H M b- H M b = hA(ojB)(cr^ n r + 
ctk},) is an effective field-mediated, dipole interaction cou¬ 
pling resonant atomic transitions. Here we assume B and 
M to be resonant through transitions at frequency 
and all their dipoles to have the same magnitude and to 
lie along the line joining the two atoms, and oriented from 
B to M. HA(uj) is the effective interaction strength and 
a (« r ) is the lowering operator of the body (of the res¬ 
onant transition of the machine). Hmb originates from 
the correlations of the fluctuations of atomic dipoles due 
to the common field. 

The derivation of the master equation Q has been 
performed under the Markovian and rotating wave ap¬ 
proximations. It involves the average photon number 
n(w,T) = 1 /( e hu / k BT _ -Q frequency u> and temper¬ 
ature T and the two functions atw(S) which encompass 
all the properties of the environment, such as the dielec¬ 
tric properties of the slab and the correlation functions 
of the field. For their explicit expressions we refer the in¬ 
terested reader to [24] . The dissipative effects due to the 
atom-field coupling are accounted for by the dissipators 
Dk with expressions 

Db(p) = T+(uj B )(vpcr i - \{<? ] o,p}) 

+ r s(ws)((T t pcr- ijcrcT^p}), (5) 

Dn{p) = F +(w n )(K„p/4 - 

+ r n k) P]) , (6) 

Dd(p) =T+(u B )[firP^ - ^{V/CtuP}) 

+ Tf(uj B )^ntpcr- ^{er/ 4 , p}) + h.c. 


( 7 ) 
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where u>d = ub- One recognises standard local dissipa¬ 
tion terms ( Db and D n ), each associated to the degrees of 
freedom of a well-identified atom, and non-local dissipa¬ 
tion ( Dd ) which describes energy exchanges at frequency 
wb of the atomic system as a whole with its OTE envi¬ 
ronment, not separable in machine or body contributions, 
its action involving degrees of freedom of both atoms in a 
symmetric way. The parameters (the rates of the 

dissipative processes of absorption and emission of pho¬ 
tons through local or non-local interactions) depend on 
local or non-local correlations of the field in the atomic 
positions, which in turn are functions of the temperatures 
Ts and TV and the dielectric properties of the slab S as 

' 0 } n = [1 + n(u,T w )]a l w (u) 

1 i 

+ [1 + n(w,r g )]ag(c<;), (8) 

=n(w,r M ,)aV(w)* +n(u;,T s )4(a;)*, (9) 

where T°(a;) = |di| 2 w 3 /(3for£oc 3 ), for i = 1,2,3, H, is 
the vacuum spontaneous emission rate of the i- th atomic 
transition having a dipole moment d 7 , and T 3 (o;) = 
(w)T°(w). Thanks to the functional dependence of 
these parameters on the frequency and on the position 
of the atom, and to the critical behaviour shown in cor¬ 
respondence to the resonance frequency us of the slab 
material, thermodynamic tasks become achievable. To 
simplify the notation, in the rest of this work the explicit 
w-dependence in all the Ts will be omitted. 

III. THERMODYNAMICS OF THE SYSTEM 

After having introduced all the dynamic effects char¬ 
acterizing the atomic system, we want in this Section 
to introduce some quantities which will characterize the 
machine tasks and functioning. 


are not simply blackbody spectra as they have 
their own transition-dependent Purcell factor [26] . 
In this framework we study thermodynamic effects of 
stationary heat fluxes between M and B , mediated and 
sustained by the OTE environment. To characterize the 
effects of these fluxes a second kind of temperature has 
to be introduced. Indeed, as much as the environmen¬ 
tal temperatures characterise the thermodynamics of the 
OTE field, we need a second parameter to describe the 
energetics of atoms. In particular, atoms exchange en¬ 
ergy under the form of heat with their surroundings by 
emitting photons through one of their transitions. This 
means that the possibility of such heat exchanges is re¬ 
lated to the distribution of population in each atomic 
level. Note that, from the very definition of T„, the en¬ 
vironmental temperature depends on how the field tends 
to distribute atomic population in each pair of levels, due 
to the presence of the ratio T+ /T~. A transition is there¬ 
fore in equilibrium with its effective local environment if 
and only if its two levels |a) and \b) are populated such 
that Pa/Pb = T+/T-. If not, the field and the atom will 
exchange heat along such a transition until such a ratio 
is reached. This suggests to introduce a second tempera¬ 
ture, hereby referred to as population temperature , which 
for a transition of frequency u n (n = 1, 2,3, B) is defined 
as 

n _ . 

" k B \n{p“/p* n y { > 

Pn ( Pn ) being the stationary population of the ground 
(excited) state of the n-th transition. The result of a 
stationary thermodynamic task on the body, be it re¬ 
frigeration, heating or population inversion, is then to 
modify its population temperature 9b- 


B. Heat fluxes 


A. Environmental and population temperatures 


In order to describe the machine thermodynamics, it is 
convenient to introduce two kind of temperatures. A first 
one characterizes the action of the field on the atoms: 
it has been shown 26] that the atom-field interaction 
can be effectively rewritten as if each atomic transition 
felt a local equilibrium environment whose temperature 
depends on the transition frequency, on the properties of 
the slab and on the slab-atom distance z. These effective 
environmental temperatures depend on the rates T„ as 


T n = T(u n ) 


hu n 

k B in(r+/r^)’ 


( 10 ) 


with n = 1,2,3, X,d. It is important to stress 

here that, despite these effective environments can 
be characterised by a temperature, their spectra 


The condition Xi = 9i is satisfied only if detailed bal¬ 
ance ( Pi/p\ = T+/T“) holds. It can be proven that 
detailed balance can be broken in a three level atom in 
OTE fields. As a consequence the machine M produces 
non-zero stationary heat fluxes with B and the field envi¬ 
ronment, one for each dissipative process D n in the ME 
©• These fluxes, following the standard approach in the 
framework of Markovian open quantum systems 128] , are 
given as Q n = Tr [H at D n p] , where p is here the station¬ 
ary atomic state and H at is a suitable atomic Hamilto¬ 
nian which can be Hm , H B or Hm + Hb depending on 
which part of the atomic system the heat flows into. Note 
that this definition implies an outgoing heat flux to be 
negative. 

Following their definition, these heat fluxes depend 
both on the field properties (through the structure of 
the dissipators D n ) and on the properties of the atoms 
through their stationary state. This dependence, for the 
local dissipators, can be put under the very clear ther- 
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modynamic form 

Qn = K n (e^ - e 5 ^) ~ C n (9 n )(T n - 9 n ), (12) 


where K n > 0, C n (9 n ) is a positive function of 9 n (and of 
other parameters such as the frequency of the transition) 
and the second approximated equality holds in the limit 
9 n ~ T n . Equation (12) shows that the direction of heat 
flowing is uniquely determined by the sign of the differ¬ 
ence T n — 9 ni matching the thermodynamic expectation 
that heat flows naturally from the hotter to the colder 
body and strengthening the physical meaning of 9 n . 

There being no time-dependence in the Hamiltonian of 
the model, the first law of thermodynamics at stationar- 
ity for the total atomic system comprises only heat terms 
and assumes the form 


without the effect of the machine, would thermalise at 
the local environmental temperature (9b = Tn)-, cor¬ 
responding to Pe/Pg = r^/r^g. This temperature is 
necessarily constrained within the range [Tw, YV] (21)]. 
Due to the particular form of the master equation pi) , in 
which all collective atomic terms involve only resonant 
atomic transitions, in the non-resonant subspace the col¬ 
lective atomic state will be diagonal in the eigenbasis of 
Hb + Hm- This is due to the fact that local dissipation 
in Eq. p]) of Section [II] induces a thermalisation with 
respect to the free atomic Hamiltonians. On the other 
hand, in the resonant atomic subspace of the eigenbasis 
of Hb + Hm spanned by the states |g), |e) of B and the 
two states |0 r ), |l r ) of the transition of M at frequency 
ujb > the most general form of the atomic stationary state 
is 


3 

Qb + Qn + Qd — 0. (13) 

n= 1 

In addition, energy is exchanged between the machine 
and the body thanks to their field-induced interaction 
Hmb- In Appendix A, following the general scheme de¬ 
veloped in [29], we show such an exchange to be un¬ 
der the form of heat. Seen by M such a flux is Q r = 
— iTv^HmIHmb, Ps])/ft while as expected B sees the flux 
— Q r . By introducing the explicit expressions for Hm 
and Hmb , one can obtain a particularly simple form for 
Qr as 

Q r = 2SwBA(w B )((7 t Hr)-, (14) 

where (cr' n r )_ = ^(<j'K r — cm' r ). In an analogous way, by 
employing Eq. |7|, one can evaluate the change in inter¬ 
nal energy of M due to the non-local heat flux exchanged 
by the atomic system with the OTE environment, given 
by Qd = Tr [H M D d p\. It is 

Qd = -fiw B Re|(cr4)[ r d - ( r d)*]}- ( 15 ) 

Finally, the change in the internal energy of B due to 
the same effect, Tt[HbDcip\ , is given by same expression 

©• 

Fig. |l(b)| shows the full scheme of such heat fluxes for a 
particular configuration of the system. The two levels of 
B will be labeled here as |<?) and |e). Despite the two-level 
assumption might seem specific, it has been shown in 
various contexts [ 3127 ] that quantum thermal machines 
only couple to some effective two-level subspaces in the 
Hilbert space of the body they are working on. A two- 
level system is therefore the fundamental building block 
of the functioning of quantum thermodynamic tasks. 


IV. COHERENCE-DRIVEN MACHINE TASKS 

The main result of this paper is the possibility to drive 
the temperature 9 B of the body outside of the range de¬ 
fined by the external reservoirs at Tw and Tg- The body, 


'Pe i r 0 0 0 \ 

0 Pe 0 r C r 0 I 

0 c* Pgi r 0 

. 0 0 0 p g0 J 


(16) 


A coherence c r is present in the decoupled basis between 
the two atomic states |gl r ) and |e0 r ) having the same 
energy. 

Note that the temperature 9b of the body increases 
monotonically with the ratio p e /p g - By tracing out the 
machine degrees of freedom from the master equation Q, 
one obtains a diagonal state with stationary populations 
p g and p e of the body B. Be now the flux Q r seen 
by B. Then the expressions for heat fluxes exchanged by 
B with its surroundings are 


Qr = —t{[Hb,HmbX), 


Qd=T d (<7 t H B K r ) - -{{Hb,^ Kr}) 
{oH b kI) - ^({# b , o -4}) 
Qb = r B {o*H b o) - )^{{H b ,o*<j}) 


( ctH b < 7 t ) - -({ H B ,aa f }) 


(17) 


+ C.C., (18) 


(19) 


where the mean values are evaluated over the stationary 
state of the total system. Exploiting its general form 
( p~6] ) , it is just a matter of straightforward calculations to 
evaluate all the mean values above. Imposing the sum of 
( p~8]) a nd ( fl9] ) to vanish (first law for B 1 analogous 
to Eq. (|13|)), one obtains 


Pe = T b - A(ojb) 

Pg Tg + A (wb) 


( 20 ) 


where 

A(w B ) = 2A(w B )Im{c r } + Re{c r [r+- (Tj)*]}. (21) 


Note now that, thanks to Eq. ( [l6| ), 44) = c r 
and (cr'K r )„ = Im(c r ), such that the first term 
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in A(wg) stems from the resonant heat — Q r = 
— 2fiwsA(wB)Im(c T .) exchanged with the machine, while 
the second is due to the non-local heat flux Qd- Eqs. 
(201 and (21) show that the thermal machine works only 
if a stationary quantum coherence c r is present. Remark¬ 
ably, it can be shown [ 30] that quantum discord m (a 
key measure of purely quantum correlations) is a mono¬ 
tonic function of the absolute value of the coherence c r in 
our system. Differently from previous studies he here 
discord between M and B is a necessary condition for any 
thermodynamic task, and represents a resource the ma¬ 
chine can use through the two different processes Q r and 
Qd- Eq. p0| means that a quantum coherence between 
machine and body modifies the stationary temperature 
of the body with respect to Tp. This modification is re¬ 
ported in Fig. [2] where the behaviour of 9 b as a function 
of the slab-atoms distance z is shown for two different 
slab thicknesses S. Four possible regimes can be singled 
out: both during refrigeration (9b < Tg) and heating 
(9b > Tb ), 9b can be either driven outside of the range 
[T\y,Ts] (strong tasks) or kept within it (light tasks). 
As a limiting case of strong heating, the body can be 
brought to infinite temperature (p e = p g ) and, further 
on, to negative ones, producing population inversion. As 
one can easily see from Fig. [2j the physics behind the 
absorption tasks is enclosed in the strong sensitivity of 
the population temperature 9b of the body to the popu¬ 
lation temperature 9m of the machine along the resonant 
transition when the body is not present. 


A. Optimal conditions for thermodynamic tasks 

It is shown in Fig. [2] that the machine has a very high 
thermal inertia, such that the body, when put into ther¬ 
mal contact with the machine having a certain tempera¬ 
ture 9Mi thermalizes with it and 9b — 9m- Fig. [3] shows 
the mechanism the machine uses to modify its population 
temperature 9m in absence of the body, thanks to the dif¬ 
ferent environmental temperatures each of its transition 
feels. This drives M out of detailed balance condition 
and allows M to keep its resonant transition tempera¬ 
ture almost constant. We label here the three transitions 
of the machine as high frequency (u>h), average frequency 
(w a ) and low frequency (u>i), one of which (suppose here 
oji, connecting states |0 r ) and |l r )) is resonant with ug- 
For simplicity, let us focus on refrigeration only, which 
we suppose to happen either through transition 2 (con¬ 
necting first and the second excited states), since in this 
configuration the high-frequency transition 3 is always 
used by an absorption refrigerator to dissipate heat into 
the environment [16]. As shown in Fig. [2j to obtain a low 
9b the resonant machine transition must be made cold. 
This is achieved by reducing the ratio pi r /po r , which in 
turn happens when: 

(a) the effective environmental temperature Th felt by 
the high frequency transition is very cold. In this way 
the environment contributes in increasing the population 


10-' 1 5 10 io 2 



FIG. 2. Stationary temperature Ob of the body (solid black 
line), machine resonant temperature 6 m in the absence of B 
(dotted pink line) and local-environment temperature Tb felt 
by the body (dot-dashed green line) versus z. The slab is 
made of sapphire and kept at Ts = 500 K while Tw = 300 K. 
The machine transition frequencies are u\ = 0.9 us, wb = 
cd 2 = 0.1 ujs and U 3 = us, us = 0.81 ■ 10 14 rads -1 being 
the first resonance frequency of sapphire (the optical data for 
the dielectric permittivity of the slab material are taken from 
[32]). The two atoms are placed at a distance r = 1 pm from 
each other. Panel (a): numerical data for a semi-infinite slab. 
Light and strong refrigeration are achieved in this configura¬ 
tion. Panel (b): same quantities for a slab of finite thickness. 
The plotted functions are in this case —1/0b, — 1 / 6 m and 
— 1 /Tb (left vertical scale), with the same color code as be¬ 
fore. The population inversion corresponds to divergent 0b 
and 6 m ■ The corresponding value of temperature can be read 
on the right vertical scale. All tasks are in this case obtained. 


of the ground state of M at the expenses of the popula¬ 
tion of its most energetic state. The resonant transition 
involves necessarily one of these two levels, and in both 
cases the effect of the high frequency transition helps re¬ 
ducing pi r /p 0r ] 

(b) the effective environmental temperature felt by the 
average transition is very hot. This, following the same 
idea, would either mean reducing the population of pi or 
increasing the one of po, thus reducing pi r /po r - 
When these two conditions are met, the machine can 
always redistribute its populations such that the ratio 
Pi r /Po r can be kept low and almost unaffected by the 
presence of another atom. The advantage of the OTE 
field configuration is that the effective field temperatures 
can be manipulated through a wide set of parameters 
involving z, S, Tw and Tg- In particular, the role of 
the resonance of the slab material is crucial [H], as ex- 
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FIG. 3. Conditions for refrigeration: effective rate temper¬ 
atures Th, T a and 7) of the local environments felt by the 
three machine transitions, for uih = ws, ui a = 0.8 tos and 
u>i = 0.2 cus (cos = 0.81- 10 14 rad s^ 1 ) versus the machine-slab 
distance z in the absence of the body. The slab and walls 
temperatures are Ts = 200 K and Tw = 300 K, and the slab 
thickness is <5 = 0.05 /im. As the plot shows, the transition 
having the same frequency as the slab resonance is much more 
strongly affected by the field emitted by the slab, such that 
its rate temperature is kept much lower than T a . This pro¬ 
duces a mechanism, shown in the inset, according to which 
excitations (yellow dots) are transferred to the intermediate 
level of the machine and removed from its upper one. This 
in turn drives the population temperature 0 i of the transition 
at uji (transition 2 in the example) to values lower than the 
one Ti of its local environment, allowing the machine to re¬ 
frigerate objects. In this configuration, introducing a body B 
at 2 = 1 /im from the slab and r = 1 /im from M, one obtains 
0 B = 160 K < T s . 


plained in the caption of Fig. [3j In the case Ts < Tw, 
transitions strongly affected by the field emitted by the 
slab feel a cold local environment. Moreover, provided 
ui a is far enough from uig, one can at the same time have 
T a — Tw- By this mechanism, M can change the tem¬ 
perature 9b, bringing it to values far outside the range 
[T S ,T W \. 


We stress here that the difference between light and 
strong tasks is a fundamental one: better than light 
tasks could in principle be done by direct connection 
of the body to one of the two real reservoirs at Ts 
or Tw, while strong tasks can not be achieved by a 
simple thermal contact with anything in the system. 
A (ujb) strongly depends on the slab-matter system dis¬ 
tance 2 and on the external temperatures through c r , 
A and r±. One can thus engineer one or many of these 
regimes at will as shown in the functioning-phase diagram 
of the machine in Fig. [4]for a fixed thickness <5 = 0.05 /im. 
All the strong and light functioning phases of the machine 
are found as a function of both Tw — Ts and 2 . 



z (/tm) 


FIG. 4. Functioning phases of the absorption machine versus 
the atoms-slab distance 2 and AT = Tw — Ts- The sapphire 
slab has a thickness 5 = 0.05 /jm and its temperature is con¬ 
tinuously changed in the range [30 K, 570 K[. The walls are at 
Tw = 300 K. Strong refrigeration (blue areas), strong heating 
(red areas), light refrigeration (cyan), light heating (orange) 
and population inversion (green) can all be obtained. 


V. EFFICIENCY AND CARNOT LIMIT 


Consider now the refrigerating regime in which the ma¬ 
chine extracts heat from the body through the transition 
2. The scheme of h eat fluxes is then exactly the one 
depicted in Fig. |l(b)[ The efficiency of this process is 


Vref — 


Qr 

Qi + Q2 


( 22 ) 


due to the fact that Q r is the power produced by the 
machine, which absorbs energy from its surroundings 
through transitions 1 and 2 (the equivalent of a work 
input) while uses transition 3 to dissipate part of the ab¬ 
sorbed energy after use (the equivalent of the spiral in a 
normal fridge). The corresponding Carnot limit rfc e f can 
be obtained by analysing the machine functioning in its 
reversible limit (zero entropy production). The instanta¬ 
neous entropy production rate r for quantum systems is 
defined as [251 


f 7=--5(p(t)||p S ‘)>0, 


(23) 


where S(p(t)\\p st ) is the so-called relative entropy [55] , 
never increasing in time under a Markovian dynamics. 
Following [16], one can apply equation (23) term by term 


to each dissipator in the master equation thanks to the 
fact that they all are under a Markovian form. One thus 
obtains 


Tr 


(D fc p st )lnpf 


> 0 , 


(24) 
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where p j? ( , k = 1, 2,3, B , d is the kernel (stationary state) 
of the single dissipator D The 3 local dissipators D n 
of the machine and the local dissipator Db of the body 
induce stationarity under the standard Gibbs form at the 
effective environmental temperature, diagonal in the free 
atomic Hamiltonian basis. The nonlocal dissipator Dd , 
in the case studied here where the dipoles of B and M 
lie along the line connecting the two atoms (and more in 
general when T± £ R), has the same kernel at environ¬ 
mental temperature Td , local in the degrees of freedom 
of M and the B. Introducing these single-dissipator sta¬ 
tionary states into equation (24) one obtains 


Q i 


Q2 

T2 


Q 3 

T :i 


Qb 

Tb 


2 Q d 

T, 


< 0, 


(25) 


which is a form of the second law at stationarity for our 


system. With the help of the first law in Eq. (131 of Sec¬ 
tion |III[ the known property of three-level atomic heat 
fluxes |T3| \Q n /Qm\ = u n /LUm Vto, n = 1,2,3 (where 
Q n is the total flux along the n-th transition) and the 
fact that in refrigeration T 3 < T 2l Td and T 2l Td < T\ 
(as commented in Section IV A) and under the condition 


Qd < 0 (other cases can be treated analogously), one ob¬ 
tains from (25) another first degree inequality. This has a 


non-trivial solution only if Td > T 2 , from which a bound 
on the efficiency in Eq. (22) can be obtained as shown 


in Appendix B. Such a bound depends only on the three 
frequencies of the machine and the temperatures of the 
effective local and non-local environments. In the case 
of refrigeration along transition 2 the Carnot efficiency 
assumes the form 

1 T 2 Td(Ti—T 3 ) . Ul 2 T 2 (T d ~T 3 ) rp rp 

2T 1 T 3 (T d -T 2 ) "T" 2u>! T 3 (T d -T 2 )’ 11 ^ ^ J 2, 

if T d < T 2 . 

(26) 



A. Efficiency at maximum power 

An important figure of merit for the realistic func¬ 
tioning of any thermal machine is how close to its 
Carnot limit it works when delivering maximum power 
(i.e., when Q r is maximised). Many bounds are known 
for different setups, limiting the efficiency at maximum 
power rf 71 to some fractions of rp [16, 34] . Remark¬ 
ably our structured OTE environment allows for refrig¬ 
eration tasks with rp 1 much closer to rp than the bound 
known for quantum absorption machines m based on 
ideal blackbody reservoirs, reading for our system rf 71 < 
0.75 r\ c . This is exemplified in Fig. [ 5 ] for a particu¬ 
lar configuration of the model. The blue triangles (left 
vertical scale) represent the ratio r\ re f /Vrlj 1 plotted ver¬ 
sus u 3 while keeping fixed uj 2 = lob = 0.1 us- The red 
dots (right vertical scale) are the power Q r plotted ver¬ 
sus the same quantity, while the red dashed line is the 
machine-body discord (right vertical scale). It is clear 
that the power is maximised at w 3 = 1.05 wg, corre¬ 
sponding to ip e f ~ 0.89 ?y ^ e f. Qr starts decreasing, as 



^3 Ms 


Qr 


Discord 


FIG. 5. The ratio r/ref /r/ref (blue triangles, left vertical axis), 
the power of the machine ( Q r , red dots, right vertical axis in 
units of 10~ 14 /rJ/s) and machine-body discord (dashed red 
line, right vertical scale in units of 10 -4 ) versus the scaled 
machine transition frequency LJs/ojg. The sapphire slab is 
semi-infinite and at Ts = 395 K, while Tw = 125 K and 
z = 4.8 pm. The transition frequency of the body is fixed 
as O.lws (o>s = 0.81 • 10 14 rads” 1 ), resonant with transi¬ 
tion 2 of the machine. The maximum power is reached for 
UJ 3 = 1.05 us, corresponding to rjre.f /rp,.f ~ 0.89. Remark¬ 
ably, discord shows a sharp resonance peak, similarly to Q r . 


classically expected when the efficiency approaches rj c , 
around UJ 3 ~ 0.9 wg, but suddenly increases again when 
0 J 3 approaches wg. This behaviour is due to the fact that, 
when one atomic transition is resonant with the charac¬ 
teristic frequency of the slab material, the atomic popu¬ 
lations are strongly affected by the field emitted by the 
slab. Hence the not-black-body nature of the total field 
become crucial (e.g., the atomic decay rate is no longer 
proportional to w 3 ), allowing to overcome bounds set by 
the blackbody physics. The role of discord as machine re¬ 
source is clearly shown here, where discord at resonance 
has a sharp peak leading to the high-power performance 
of M. 

One could wonder whether such an exceptionally high 
efficiency at maximum power is only seldomly attained 
for the kind of machines described here. To answer 
such a question on quantitative bases, we performed 
a random sampling of over 2 • 10 4 thermal machines, 
all delivering thermodynamic tasks on the same fixed 
body. In the simulations performed and reported in 
Fig. § the machines work as a quantum refrigera¬ 
tor delivering strong refrigeration using a semi-infinite 
slab. In this sampling, the machine-slab distance z has 
been, for each machine, randomly drawn in the range 
[0.9 /rnr, 100 /irn] , the walls temperature has been selected 
randomly in Tw £ [50 K, 500 K] and, for each value of 
Tw , the slab temperature has been chosen at random 
in Ts £ [Tw,Tw + 500 K]. The internal structure of 
the body is kept fixed during the simulations, with a 
frequency lob = 0.1 u>s resonant with the transition 2 
of M. For each machine thus generated, we have then 
maximised the delivered power by modifying the two 

















FIG. 6. Statistical occurrence of ratios rf 1 /rf for a random 
sampling of 2-10 4 thermal machines, always in resonance with 
the same body B. For each machine, z has been randomly 
generated in the range [0.9/xm, 100/xm], Tw £ [50K,500K] 
and, for each value of Tw , Ts £ [TV, TV + 500 K]. The inter¬ 
nal structure of the body is kept fixed during the simulations, 
with a frequency ojb = 0.1 ws resonant with the transition 
2 of M. The maximisation is performed over every possible 
value of toi £ ( 0 J 2 ,tos) and L 03 compatible with the condition 
oji + u >2 = UJ 3 . Around 50% of machines thus generated have 
rj m > 0.75 rf. 


other machine frequencies over every possible value of 
to 1 £ (u> 2 ,ujs) and 0 J 3 compatible with the condition 
wi + u >2 = 0 J 3 . Finally, once obtained the configuration 
corresponding to the maximum power, we have computed 
the efficiency of the process. Fig. [6] shows the histogram 
of the distribution of the ratio rf 1 /rf of efficiency at 
maximum power to the corresponding Carnot efficiency 
in the interval [0,1] within these 2 • 10 4 random refrigera¬ 
tors. It is remarkable that around 50% of these machines 
work at maximum power with efficiencies higher than the 
bound 0.75 rf in [IB] and that none of them have been 
found to work at maximum power with efficiencies lower 
than 0.6r] C . Moreover, as can be clearly seen in Fig. 
[o] a small but non-negligible fraction of them can reach 
rf ~ 0.98 rf. 


VI. CONCLUSIONS 

This work introduces a new realization of a quantum 
thermal machine using atoms interacting with single non- 
equilibrium electromagnetic fields. By simply connect¬ 
ing two thermal reservoirs to macroscopic objects , their 
radiated held allows the atomic machine to achieve all 
quantum thermodynamic effects (heating, cooling, pop¬ 
ulation inversion), without any direct external manipula¬ 
tion of atomic interactions. This overcomes the usual dif¬ 
ficulty of connecting single transitions to thermal reser¬ 
voirs, in a realistic and simple configuration where the 
field-mediated atomic interaction modifies at will station¬ 
ary inter-atomic energy fluxes. 


Despite the environmental dissipative effects, atoms 
share steady quantum correlations [24[ 1.851 which we 
showed to be necessary for one atom to deliver a ther¬ 
modynamic task on the other, uncovering genuinely non- 
classical machine functioning. These particular features 
affect the tasks efficiency, which can be remarkably high 
also at maximum power, defying the known bounds for 
quantum machines based on ideal and independent black- 
body reservoirs thanks to the fundamental effect of the 
resonance with the real material of which the slab is 
made. Moreover, such a remarkably high efficiency at 
maximum power is strongly connected to the presence of 
a peak in quantum correlations between the machine and 
the body, which represent the resource the machine uses 
for its tasks. 

These results tackle major open problems on quantum 
thermal machines, paving the way for an efficient 
quantum energy management based on the potentialities 
of non-equilibrium and quantum features in atomic-scale 
thermodynamics. 
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Appendix A: Resonant heat flux 


In this appendix we demonstrate that the resonant 
energy exchange between M and B due to the field- 
mediated coherent interaction Hmb consists only of heat. 
Following the approach of [29], the dynamics of the sole 
M induced by the Hamiltonian interaction Hmb com¬ 
prises in general an Hamiltonian and a dissipative part 
and can be written as 


Pm — ~\Hm + H°m , pm\ + Dmb(p), (Al) 


where Dmb is a non-unitary dissipative term for M due 
to the interaction with B , which depends however on the 
total state p = pmb because, in general, the two subparts 
are correlated. H^f is a renormalised free Hamiltonian 
of subsystem M due to the interaction with B. Defining 
the two marginals Pm(B) = Ti 'b(m)P and the correlation 
operator Cmb = P — Pm ® Pb, it is shown in [29] that 

Hm = Ti\b (h mb ()-m ® Pb )) > (A2) 

Dmb {p) = -*Tr B Cmb]^ ■ (A3) 


Introducing H^i as the part of Hm which commutes 
with Hm and Hm 2 w hi c h does not, directly from equa¬ 
tion (Al) one has, for the internal energy of M Um = 
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Tr((H M + H$)p), 

Um = Ttm ({Hm + )Dmb{p)^ + Ti'm(^?miPm) 

- *Tr M ([% + Hlf ,, p M ) • 

(A4) 

It is custom to identify heat terms as the ones producing 
a change in the entropy of a subsystem: all the rest is 
identified as work W. Eq. (A4) can then be split in 


Qm = Tr M (( Hm + H^DMBip )), (A5) 

Wm = Tr M (h$ iPm ~ i \Hm + fflm ^A/2] Pm)(- A6) 

Introducing the symbols c ^ = (i|/?M|j) (* 7^ j) for the 
coherences of the marginal p S M (different then from the 
coherence c r introduced in Eq. (16) of Section |IV| which 


is a two-atom coherence), equation (A2) becomes 

a Re(c“). 


(A7) 


By tracing out the machine or the body degrees of free¬ 
dom from equation (16), one can prove that the two sta¬ 
tionary marginals p ^ and pff are always diagonal in the 
eigenbases of their respective free Hamiltonians, so that 
c m = 0- No renormalisation to the machine Hamilto¬ 
nian comes therefore from the interaction with B , which 
means that equation ( |A6| ) vanishes, proving that no work 
is involved in machine-body energy exchanges. As for the 
heat, considering that [Hmb , Pm ® Pb\ — 0, Eq. (A5) re¬ 
duces to 


Qm = — iTr M (^H M Th:B[H M B, p st ]^ = Qr 


(AS) 


with the same Q r given in Eq. (14). 


Appendix B: Carnot limit 


In this appendix we deduce Eq. (26) of Section [V] for 
the Carnot efficiency in refrigeration along transition 2, 
and under the condition Q d < 0. In addition to Eqs. 
(13) and ([25]) , the condition \Q n /Qm\ = w n /w m gives for 


n = 1 and m = 2 the following 

Qi 


Q 2 + Qr + Qd 


Wl 

U >2 


(Bl) 


Solving Eqs. (13) and (Bl) for Q 3 and Q d and using 


these solutions into (25) one obtains for Q 


Q r if Q 1 


T d T 2 

T d -T 2 


h(i+^ 


T. 


Wl 


) 1 

UJ 2 1 

' Ti 

wi T d 


~Q 2 (B2) 


which, used in Eq. (22) of Section [VJ gives a bound 
on r/ re f as a function of Q\ and Q 2 . Finally, using the 
fact that such a bound is a decreasing function of Q 2 , 
one obtains the Carnot efficiency as the limit for Q 2 —► 
0, which turns out to be independent on Q\ and gives 
ultimately the first line of Eq. (26). On the other hand, 


in the case T 2 < T d , one can not obtain anything like Eq. 
(B2) and the only possibility for the machine to work 


without producing entropy is therefore to have vanishing 
heat flux from/to the body. This means Q 2 = Q d = 0 
which, inserted in the expression for the efficiency and 
using again \Q n /Qm\ = ea n /uj m leads to the second line 
of Eq. ([26b. 
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